WebCM: A Web-Based Platform for Multiuser Individual-Based Modeling of Multicellular Microbial Populations and Communities.

WebCM is a web platform that enables users to create, edit, run, and view individual-based simulations of multicellular microbial populations and communities on a remote compute server. WebCM builds upon the simulation software CellModeller in the back end and provides users with a web-browser-based modeling interface including model editing, execution, and playback. Multiple users can run and manage multiple simulations simultaneously, sharing the host hardware. Since it is based on CellModeller, it can utilize both GPU and CPU parallelization. The user interface provides real-time interactive 3D graphical representations for inspection of simulations at all time points, and the results can be downloaded for detailed offline analysis. It can be run on cloud computing services or on a local server, allowing collaboration within and between laboratories.

Flapjack: Data Management and Analysis for Genetic Circuit Characterization.

Flapjack presents a valuable solution for addressing challenges in the Design, Build, Test, Learn (DBTL) cycle of engineering synthetic genetic circuits. This platform provides a comprehensive suite of features for managing, analyzing, and visualizing kinetic gene expression data and associated metadata. By utilizing the Flapjack platform, researchers can effectively integrate the test phase with the build and learn phases, facilitating the characterization and optimization of genetic circuits. With its user-friendly interface and compatibility with external software, the Flapjack platform offers a practical tool for advancing synthetic biology research.This chapter provides an overview of the data model employed in Flapjack and its hierarchical structure, which aligns with the typical steps involved in conducting experiments and facilitating intuitive data management for users. Additionally, this chapter offers a detailed description of the user interface, guiding readers through accessing Flapjack, navigating its sections, performing essential tasks such as uploading data and creating plots, and accessing the platform through the pyFlapjack Python package.

Genetic Network Design Automation with LOICA.

Genetic design automation (GDA) is the use of computer-aided design (CAD) in designing genetic networks. GDA tools are necessary to create more complex synthetic genetic networks in a high-throughput fashion. At the core of these tools is the abstraction of a hierarchy of standardized components. The components' input, output, and interactions must be captured and parametrized from relevant experimental data. Simulations of genetic networks should use those parameters and include the experimental context to be compared with the experimental results.This chapter introduces Logical Operators for Integrated Cell Algorithms (LOICA), a Python package used for designing, modeling, and characterizing genetic networks using a simple object-oriented design abstraction. LOICA represents different biological and experimental components as classes that interact to generate models. These models can be parametrized by direct connection to the Flapjack experimental data management platform to characterize abstracted components with experimental data. The models can be simulated using stochastic simulation algorithms or ordinary differential equations with varying noise levels. The simulated data can be managed and published using Flapjack alongside experimental data for comparison. LOICA genetic network designs can be represented as graphs and plotted as networks for visual inspection and serialized as Python objects or in the Synthetic Biology Open Language (SBOL) format for sharing and use in other designs.

A computational framework for testing hypotheses of the minimal mechanical requirements for cell aggregation using early annual killifish embryogenesis as a model.

Introduction: Deciphering the biological and physical requirements for the outset of multicellularity is limited to few experimental models. The early embryonic development of annual killifish represents an almost unique opportunity to investigate de novo cellular aggregation in a vertebrate model. As an adaptation to seasonal drought, annual killifish employs a unique developmental pattern in which embryogenesis occurs only after undifferentiated embryonic cells have completed epiboly and dispersed in low density on the egg surface. Therefore, the first stage of embryogenesis requires the congregation of embryonic cells at one pole of the egg to form a single aggregate that later gives rise to the embryo proper. This unique process presents an opportunity to dissect the self-organizing principles involved in early organization of embryonic stem cells. Indeed, the physical and biological processes required to form the aggregate of embryonic cells are currently unknown. Methods: Here, we developed an in silico, agent-based biophysical model that allows testing how cell-specific and environmental properties could determine the aggregation dynamics of early Killifish embryogenesis. In a forward engineering approach, we then proceeded to test two hypotheses for cell aggregation (cell-autonomous and a simple taxis model) as a proof of concept of modeling feasibility. In a first approach (cell autonomous system), we considered how intrinsic biophysical properties of the cells such as motility, polarity, density, and the interplay between cell adhesion and contact inhibition of locomotion drive cell aggregation into self-organized clusters. Second, we included guidance of cell migration through a simple taxis mechanism to resemble the activity of an organizing center found in several developmental models. Results: Our numerical simulations showed that random migration combined with low cell-cell adhesion is sufficient to maintain cells in dispersion and that aggregation can indeed arise spontaneously under a limited set of conditions, but, without environmental guidance, the dynamics and resulting structures do not recapitulate in vivo observations. Discussion: Thus, an environmental guidance cue seems to be required for correct execution of early aggregation in early killifish development. However, the nature of this cue (e.g., chemical or mechanical) can only be determined experimentally. Our model provides a predictive tool that could be used to better characterize the process and, importantly, to design informed experimental strategies.

Experimental Data Connector (XDC): Integrating the Capture of Experimental Data and Metadata Using Standard Formats and Digital Repositories.

Accelerating the development of synthetic biology applications requires reproducible experimental findings. Different standards and repositories exist to exchange experimental data and metadata. However, the associated software tools often do not support a uniform data capture, encoding, and exchange of information. A connection between digital repositories is required to prevent siloing and loss of information. To this end, we developed the Experimental Data Connector (XDC). It captures experimental data and related metadata by encoding it in standard formats and storing the converted data in digital repositories. Experimental data is then uploaded to Flapjack and the metadata to SynBioHub in a consistent manner linking these repositories. This produces complete connected experimental data sets that are exchangeable. The information is captured using a single template Excel Workbook, which can be integrated into existing experimental workflow automation processes and semiautomated capture of results.

Accurate characterization of dynamic microbial gene expression and growth rate profiles.

Genetic circuits are subject to variability due to cellular and compositional contexts. Cells face changing internal states and environments, the cellular context, to which they sense and respond by changing their gene expression and growth rates. Furthermore, each gene in a genetic circuit operates in a compositional context of genes which may interact with each other and the host cell in complex ways. The context of genetic circuits can, therefore, change gene expression and growth rates, and measuring their dynamics is essential to understanding natural and synthetic regulatory networks that give rise to functional phenotypes. However, reconstruction of microbial gene expression and growth rate profiles from typical noisy measurements of cell populations is difficult due to the effects of noise at low cell densities among other factors. We present here a method for the estimation of dynamic microbial gene expression rates and growth rates from noisy measurement data. Compared to the current state-of-the-art, our method significantly reduced the mean squared error of reconstructions from simulated data of growth and gene expression rates, improving the estimation of timing and magnitude of relevant shapes of profiles. We applied our method to characterize a triple-reporter plasmid library combining multiple transcription units in different compositional and cellular contexts in Escherichia coli. Our analysis reveals cellular and compositional context effects on microbial growth and gene expression rate dynamics and suggests a method for the dynamic ratiometric characterization of constitutive promoters relative to an in vivo reference.

LOICA: Integrating Models with Data for Genetic Network Design Automation.

Genetic design automation tools are necessary to expand the scale and complexity of possible synthetic genetic networks. These tools are enabled by abstraction of a hierarchy of standardized components and devices. Abstracted elements must be parametrized from data derived from relevant experiments, and these experiments must be related to the part composition of the abstract components. Here we present Logical Operators for Integrated Cell Algorithms (LOICA), a Python package for designing, modeling, and characterizing genetic networks based on a simple object-oriented design abstraction. LOICA uses classes to represent different biological and experimental components, which generate models through their interactions. These models can be parametrized by direct connection to data contained in Flapjack so that abstracted components of designs can characterize themselves. Models can be simulated using continuous or stochastic methods and the data published and managed using Flapjack. LOICA also outputs SBOL3 descriptions and generates graph representations of genetic network designs.

Flapjack: Data Management and Analysis for Genetic Circuit Characterization.

Characterization is fundamental to the design, build, test, learn (DBTL) cycle for engineering synthetic genetic circuits. Components must be described in such a way as to account for their behavior in a range of contexts. Measurements and associated metadata, including part composition, constitute the test phase of the DBTL cycle. These data may consist of measurements of thousands of circuits, measured in hundreds of conditions, in multiple assays potentially performed in different laboratories and using different techniques. In order to inform the learn phase this large volume of data must be filtered, collated, and analyzed. Characterization consists of using this data to parametrize models of component function in different contexts, and combining them to predict behaviors of novel circuits. Tools to store, organize, share, and analyze large volumes of measurement and metadata are therefore essential to linking the test phase to the build and learn phases, closing the loop of the DBTL cycle. Here we present such a system, implemented as a web app with a backend data registry and analysis engine. An interactive frontend provides powerful querying, plotting, and analysis tools, and we provide a REST API and Python package for full integration with external build and learn software. All measurements are associated with circuit part composition via SBOL (Synthetic Biology Open Language). We demonstrate our tool by characterizing a range of genetic components and circuits according to composition and context.

Novel Tunable Spatio-Temporal Patterns From a Simple Genetic Oscillator Circuit.

Multicellularity, the coordinated collective behavior of cell populations, gives rise to the emergence of self-organized phenomena at many different spatio-temporal scales. At the genetic scale, oscillators are ubiquitous in regulation of multicellular systems, including during their development and regeneration. Synthetic biologists have successfully created simple synthetic genetic circuits that produce oscillations in single cells. Studying and engineering synthetic oscillators in a multicellular chassis can therefore give us valuable insights into how simple genetic circuits can encode complex multicellular behaviors at different scales. Here we develop a study of the coupling between the repressilator synthetic genetic ring oscillator and constraints on cell growth in colonies. We show in silico how mechanical constraints generate characteristic patterns of growth rate inhomogeneity in growing cell colonies. Next, we develop a simple one-dimensional model which predicts that coupling the repressilator to this pattern of growth rate via protein dilution generates traveling waves of gene expression. We show that the dynamics of these spatio-temporal patterns are determined by two parameters; the protein degradation and maximum expression rates of the repressors. We derive simple relations between these parameters and the key characteristics of the traveling wave patterns: firstly, wave speed is determined by protein degradation and secondly, wavelength is determined by maximum gene expression rate. Our analytical predictions and numerical results were in close quantitative agreement with detailed individual based simulations of growing cell colonies. Confirming published experimental results we also found that static ring patterns occur when protein stability is high. Our results show that this pattern can be induced simply by growth rate dilution and does not require transition to stationary phase as previously suggested. Our method generalizes easily to other genetic circuit architectures thus providing a framework for multi-scale rational design of spatio-temporal patterns from genetic circuits. We use this method to generate testable predictions for the synthetic biology design-build-test-learn cycle.

An equilibrium model for ribosome competition.

The number of ribosomes in a cell is considered as limiting, and gene expression is thus largely determined by their cellular concentration. In this work we develop a toy model to study the trade-off between the ribosomal supply and the demand of the translation machinery, dictated by the composition of the transcript pool. Our equilibrium framework is useful to highlight qualitative behaviours and new means of gene expression regulation determined by the fine balance of this trade-off. We also speculate on the possible impact of these mechanisms on cellular physiology.

Artificial Symmetry-Breaking for Morphogenetic Engineering Bacterial Colonies.

Morphogenetic engineering is an emerging field that explores the design and implementation of self-organized patterns, morphologies, and architectures in systems composed of multiple agents such as cells and swarm robots. Synthetic biology, on the other hand, aims to develop tools and formalisms that increase reproducibility, tractability, and efficiency in the engineering of biological systems. We seek to apply synthetic biology approaches to the engineering of morphologies in multicellular systems. Here, we describe the engineering of two mechanisms, symmetry-breaking and domain-specific cell regulation, as elementary functions for the prototyping of morphogenetic instructions in bacterial colonies. The former represents an artificial patterning mechanism based on plasmid segregation while the latter plays the role of artificial cell differentiation by spatial colocalization of ubiquitous and segregated components. This separation of patterning from actuation facilitates the design-build-test-improve engineering cycle. We created computational modules for CellModeller representing these basic functions and used it to guide the design process and explore the design space in silico. We applied these tools to encode spatially structured functions such as metabolic complementation, RNAPT7 gene expression, and CRISPRi/Cas9 regulation. Finally, as a proof of concept, we used CRISPRi/Cas technology to regulate cell growth by controlling methionine synthesis. These mechanisms start from single cells enabling the study of morphogenetic principles and the engineering of novel population scale structures from the bottom up.

Orthogonal intercellular signaling for programmed spatial behavior.

Bidirectional intercellular signaling is an essential feature of multicellular organisms, and the engineering of complex biological systems will require multiple pathways for intercellular signaling with minimal crosstalk. Natural quorum-sensing systems provide components for cell communication, but their use is often constrained by signal crosstalk. We have established new orthogonal systems for cell-cell communication using acyl homoserine lactone signaling systems. Quantitative measurements in contexts of differing receiver protein expression allowed us to separate different types of crosstalk between 3-oxo-C6- and 3-oxo-C12-homoserine lactones, cognate receiver proteins, and DNA promoters. Mutating promoter sequences minimized interactions with heterologous receiver proteins. We used experimental data to parameterize a computational model for signal crosstalk and to estimate the effect of receiver protein levels on signal crosstalk. We used this model to predict optimal expression levels for receiver proteins, to create an effective two-channel cell communication device. Establishment of a novel spatial assay allowed measurement of interactions between geometrically constrained cell populations via these diffusible signals. We built relay devices capable of long-range signal propagation mediated by cycles of signal induction, communication and response by discrete cell populations. This work demonstrates the ability to systematically reduce crosstalk within intercellular signaling systems and to use these systems to engineer complex spatiotemporal patterning in cell populations.

Characterization of Intrinsic Properties of Promoters.

Accurate characterization of promoter behavior is essential for the rational design of functional synthetic transcription networks such as logic gates and oscillators. However, transcription rates observed from promoters can vary significantly depending on the growth rate of host cells and the experimental and genetic contexts of the measurement. Furthermore, in vivo measurement methods must accommodate variation in translation, protein folding, and maturation rates of reporter proteins, as well as metabolic load. The external factors affecting transcription activity may be considered to be extrinsic, and the goal of characterization should be to obtain quantitative measures of the intrinsic characteristics of promoters. We have developed a promoter characterization method that is based on a mathematical model for cell growth and reporter gene expression and exploits multiple in vivo measurements to compensate for variation due to extrinsic factors. First, we used optical density and fluorescent reporter gene measurements to account for the effect of differing cell growth rates. Second, we compared the output of reporter genes to that of a control promoter using concurrent dual-channel fluorescence measurements. This allowed us to derive a quantitative promoter characteristic (ρ) that provides a robust measure of the intrinsic properties of a promoter, relative to the control. We imposed different extrinsic factors on growing cells, altering carbon source and adding bacteriostatic agents, and demonstrated that the use of ρ values reduced the fraction of variance due to extrinsic factors from 78% to less than 4%. This is a simple and reliable method to quantitatively describe promoter properties.

Cell polarity-driven instability generates self-organized, fractal patterning of cell layers.

As a model system to study physical interactions in multicellular systems, we used layers of Escherichia coli cells, which exhibit little or no intrinsic coordination of growth. This system effectively isolates the effects of cell shape, growth, and division on spatial self-organization. Tracking the development of fluorescence-labeled cellular domains, we observed the emergence of striking fractal patterns with jagged, self-similar shapes. We then used a large-scale, cellular biophysical model to show that local instabilities due to polar cell-shape, repeatedly propagated by uniaxial growth and division, are responsible for generating this fractal geometry. Confirming this result, a mutant of E. coli with spherical shape forms smooth, nonfractal cellular domains. These results demonstrate that even populations of relatively simple bacterial cells can possess emergent properties due to purely physical interactions. Therefore, accurate physico-genetic models of cell growth will be essential for the design and understanding of genetically programmed multicellular systems.

Synthetic Biology: opportunities for Chilean bioindustry and education.

In an age of pressing challenges for sustainable production of energy and food, the new field of Synthetic Biology has emerged as a promising approach to engineer biological systems. Synthetic Biology is formulating the design principles to engineer affordable, scalable, predictable and robust functions in biological systems. In addition to efficient transfer of evolved traits from one organism to another, Synthetic Biology offers a new and radical approach to bottom-up engineering of sensors, actuators, dynamical controllers and the biological chassis they are embedded in. Because it abstracts much of the mechanistic details underlying biological component behavior, Synthetic Biology methods and resources can be readily used by interdisciplinary teams to tackle complex problems. In addition, the advent of robust new methods for the assembly of large genetic circuits enables teaching Biology and Bioengineering in a learning-by-making fashion for diverse backgrounds at the graduate, undergraduate and high school levels. Synthetic Biology offers unique opportunities to empower interdisciplinary training, research and industrial development in Chile for a technology that promises a significant role in this century's economy.

Synthetic Biology: opportunities for Chilean bioindustry and education

In an age of pressing challenges for sustainable production of energy and food, the new field of Synthetic Biology has emerged as a promising approach to engineer biological systems. Synthetic Biology is formulating the design principles to engineer affordable, scalable, predictable and robust functions in biological systems. In addition to efficient transfer of evolved traits from one organism to another, Synthetic Biology offers a new and radical approach to bottom-up engineering of sensors, actuators, dynamical controllers and the biological chassis they are embedded in. Because it abstracts much of the mechanistic details underlying biological component behavior, Synthetic Biology methods and resources can be readily used by interdisciplinary teams to tackle complex problems. In addition, the advent of robust new methods for the assembly of large genetic circuits enables teaching Biology and Bioengineering in a learning-by-making fashion for diverse backgrounds at the graduate, undergraduate and high school levels. Synthetic Biology offers unique opportunities to empower interdisciplinary training, research and industrial development in Chile for a technology that promises a significant role in this century's economy.

Computational modeling of synthetic microbial biofilms.

Microbial biofilms are complex, self-organized communities of bacteria, which employ physiological cooperation and spatial organization to increase both their metabolic efficiency and their resistance to changes in their local environment. These properties make biofilms an attractive target for engineering, particularly for the production of chemicals such as pharmaceutical ingredients or biofuels, with the potential to significantly improve yields and lower maintenance costs. Biofilms are also a major cause of persistent infection, and a better understanding of their organization could lead to new strategies for their disruption. Despite this potential, the design of synthetic biofilms remains a major challenge, due to the complex interplay between transcriptional regulation, intercellular signaling, and cell biophysics. Computational modeling could help to address this challenge by predicting the behavior of synthetic biofilms prior to their construction; however, multiscale modeling has so far not been achieved for realistic cell numbers. This paper presents a computational method for modeling synthetic microbial biofilms, which combines three-dimensional biophysical models of individual cells with models of genetic regulation and intercellular signaling. The method is implemented as a software tool (CellModeller), which uses parallel Graphics Processing Unit architectures to scale to more than 30,000 cells, typical of a 100 µm diameter colony, in 30 min of computation time.

Fast image reconstruction in fluorescence optical tomography using data compression.

We present a method for fast reconstruction in fluorescence optical tomography with very large data sets. In recent reports, CCD cameras at multiple positions have been used to collect optical measurements, producing more than 10(7) data samples. This makes storage of the full system Jacobian infeasible, and so data are usually subsampled. The method reported here allows use of the full data set, via image compression methods, and explicit construction of the (small) Jacobian, meaning optimal inversion methods can be applied, and thus leading to very fast reconstruction.